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I. INTRODUCTION 



The main cosmological constant problem (CCPl) can be phrased as follows (see, e.g., the 
review why do the quantum fields in the vacuum not produce naturally a large absolute 
value for the cosmological constant A in the Einstein field equations or, practically, why is 
the measured value of \A\^^'^ very much smaller than the known energy scales of elementary 
particle physics ? One possible solution relies on the so-called g-theory approach [^44] , which 
provides a compensation- type solution of CCPl by self-adjustment of the q{x) field. This 
q field, considered to describe ultrahigh-energy microscopic degrees of freedom, must be of 
a very special type, being relativistic and conserved in the equilibrium state (Minkowski 
spacetime). In fact, the equilibrium value go is constant over spacetime. This property 
allows for the study of the macroscopic equations in terms of q^, which, in particular, give 
a vanishing gravitating vacuum energy density (provided there is no external pressure): 
pviqo) = Ao = 0. 

In this way, the zero-point energies of the standard-model fields in the equilibrium vacuum 
state can be compensated exactly by contributions from the microscopic degrees of freedom 
at a higher energy scale, without need to know the detailed microscopic theory. Having 
provided a possible explanation of the vanishing cosmological constant Aq in the equilibrium 
state, the next task is to explain the measured small but nonzero value of the effective 
cosmological constant A^s in the present expanding (nonequilibrium) Universe. The search 
for the explanation of this last number, Aefr, has been called the second cosmological constant 
problem (CCP2), even though the term 'puzzle' is perhaps more appropriate. 

In the early phase of the history of the Universe (close to the Planck epoch), the g- 
theory dynamical equations js, 0] show that the gravitating vacuum energy density pv[<l{t)] 
is rapidly relaxed to zero. What happens next depends on the details of the particle-physics 
theory, in particular, the theory at the TeV energy scale jsj. If there exist new ultramassive 
TeV-scale particles with electroweak interactions, the presence of these ultramassive parti- 
cles affects the expansion rate of the Universe, and this change of the expansion rate "kicks" 
pv{t) away from zero js, 3]- The maximal value of pv{t) is of the order of (meV)^, consistent 
with the value suggested by observational cosmology. The problem, however, is that pv(t) 
can be expected to drop to zero again if^he ultramassive particles ultimately disappear. 

Possible quantum-dissipative effects [6] may lead to a freezing of the gravitating vacuum 
energy density pv, but this does not guarantee an asymptotic approach to a standard de- 
Sitter universe (consistent with the ACDM model of the present Universe). In fact, there 
appears to be a potential mismatch [2] between the g-theory dynamical equations and those 
of standard general relativity with a nonzero cosmological constant. It is possible to modify 
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the relevant g-theory dynamical equation by hand [tI], but such a procedure is unsatisfactory. 

In this article, it is shown that it is possible to remain entirely within the framework of 
g-theory by allowing for a nontrivial interaction between, on the one hand, the g-field and, 
on the other hand, the matter and gravitational fields. The model is remarkably simple 
and has one crucial ingredient, which, ultimately, needs to be derived from the underlying 
microscopic theory (assuming that the model is relevant). The purpose of this article is to 
present the simple model and to perform a numerical calculation in order to make sure that 
there is indeed a nonzero remnant vacuum energy density, leaving an extensive discussion 
of the systematics to a future publication. 

In order to place the present work in context, the reader is referred to a recent review 
article on g-theory and the evolution of vacuum energy density in cosmology ^]. At this 
moment, it may also be helpful to comment on how the results of the present article compare 
with those of so-called quintessence models from a dynamic scalar field ay see, e.g., 
Sec. 8 of the review pjj] for further discussion and references). There are, at least, two basic 
differences. 

First, the standard quintessence models have a fundamental canonical scalar field 
0(x) and no natural explanation of CCPl (see, e.g.. Sec. VI of Ref. jl|), whereas the 
(pseudo-) scalar q{x) field is nonfundamental and special, in order to provide a possible 
solution of CCPl as explained in Ref. j4|. (Remark that quintessence models may still 
provided valuable insights into CCP2, if CCPl can be assumed to be solved.) Second, the 
"dark energy" from dynamic scalar fields has generally an equation-of-state (EOS) param- 
eter w = P/p 7^ —1, whereas the vacuum energy density of g-theory in its simplest form 
has w = —1 exactly. As regards g-theory, both points are clarified by considering explicit 
realizations, see Refs. and Sec. Ill Al 

In short, the ultimate goal is to find an explanation of both CCPl and CCP2. This article 
tries to make a modest contribution towards reaching that goal. 



II. FIELD-THEORETIC MODEL 
A. General properties 

Consider two real scalars: an ultramassive scalar field a with mass M = E^^ ~ TeV and 
a strictly massless scalar field ip. Let a now exists in A^^i identical copies and ip in N2 copies. 
Then, the scalars ipc, for c = 1, . . . , A^2 and N2 ~ 10^, may correspond to the particles of 
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the Standard ModeP and the scalars ab, for b = 1, . . . , A^i, to the particles of new TeV- 
scale physics (perhaps with Ni = N2 ^ 10^ resulting from broken supersymmetry jl2l . [fsj]). 
From now on, the indices b, c will be kept implicit by using an inner-product notation with 

a ■ a = Y,b (^bCTb and ip ■ ip = Y,ci'c '4'c- 

This article employs the framework of g-theory, possibly viewed as an effective theory 
incorporating quantum effects of vacuum- matter interactions. Following Refs. js, ^, the g- 
theory is considered to be realized via a three- form gauge field A IJ, ll5| . The macroscopic 
effective action of the relevant "ultraviolet" fields (here, A) and the "infrared" fields (here, 
g, a, and ip) is taken to be of the following form: 

S,f,,T[A, g, a, tP] = [ d^x V-detig) fc(g) R[g] + ey (g) + 4ff t[^> 9]) , (2.1a) 

q ^ e^^-'^VaA^.s/v^^d^), (2.1b) 

where R[g] is the Ricci curvature scalar obtained from the metric gap, 60/37,5 the Levi- 
Civita tensor density, and the standard covariant derivative. The effective gravitational- 
coupling parameter Kt in (12. lap is an even function of the vacuum variable q = q[A, g] and 
an indirect function of the matter fields manifesting itself as a dependence on the temper- 
ature T (see below). The energy density ev{q) is assumed to be a generic even function of 
the vacuum variable q, that is, a function different from the quadratic | corresponding to 
a Maxwell-type theory Q, 13]. Here, and in the following, set H = c = k = 1 and take the 
metric signature (—1, 1, 1, 1). 

The g-dependence of K, first introduced in Ref. j^, allows for a continuous time depen- 
dence of g in a cosmo logical context. This behavior differs from the stepwise evolution of g 



in the Brown-Teitelboim mechanism which operates via quantum tunneling |16l . Il7l |. The 
physical motivation for having K{q) is that g is considered to be one of the variables which 
characterize the quantum vacuum and that, therefore, g can be expected to determine all 
"constants" of the low-energy theory, one of which controls the gravitational coupling. The 
possible g-dependence of the "constants" in the matter Lagrange density of (12. lap is 
neglected for simplicity. As it stands, the action (12. lap corresponds to a type of Brans-Dicke 



model [18| with a special nonfundamental pseudoscalar field g. 



The nonfundamental pseudoscalar field q{x) from fl2.1bp has been called the "quinta- 



essentia" field in Ref. in order to distinguish it from the fundamental scalar field (j){x) 



^ In principle, it is no problem to extend the theory of ipc scalars to the one of the standard model with gauge 
and Higgs bosonic fields and quark and lepton fermionic fields. Actually, the fermions give the largest 
contribution to the total number of degrees of freedom in the standard model, hence the suggestive 
notation ipc for the corresponding scalars of the simplified model. 
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of typical quintessence models jsUl- The explicit realization of q via (12. ip clarifies the 
two points mentioned in the penultimate paragraph of Sec. [H q being nonfundamental and 
giving rise to a dynamical vacuum energy density with EOS parameter wy = —1- Note that, 
for this particular realization of g-theory, the mass dimension of q equals 2. For different 
realizations, q may have different mass dimensions and intrinsic parities. 



B. Specific properties 

As explained in the Introduction, the main focus is on the freezing mechanism of the 
vacuum energy density and the scalar Lagrange density is kept as simple as possible, only 
asking that it produces massive a scalars at an appropriate epoch in the history of the 
Universe (temperatures of order -Eew)- Specifically, take the following Lagrange density with 
a single quartic coupling term between the two types of scalars: 

4?r[^, = \ do.^ ■ + ^ d^a ■ d'^o ^\ {a ■ a) + (^ ■ ^P) {(t ■ a) , (2.2a) 

M = E,^, (2.2b) 

fo for T>Te,3, 

9t = I / X (2.2c) 

[go(l-{T/T,,gf) hTT<T,,g, 

where Tc^g is a critical temperature of order E^^, above which the scalar interactions are 
suppressed. The nontrivial temperature behavior of (I2.2cl) may effectively result from an 
interaction term {ip - ip) {a ■ a) j (-Eew)^ in an extended theory where a single x scalar picks 
up an expectation value at temperatures below a second-order continuous phase transition 
(see, e.g.. Sec. 4.4 in Ref. [isl). But, most likely, this phase-transition explanation of the 
postulated behavior (I2.2cp does not need to be taken literally: gT may very well have an 
entirely different origin, provided the model is relevant at all.^ 

The gravitating vacuum energy density near equilibrium (g = go) is taken to be 
quadratic 

m 

1 2 

PF(g) = e(g) -/iog= - (g-go) , (2.2d) 



^ The quartic interaction term in p.2a|) leads to radiative corrections for the low-energy theory of the 
scalars ■(/'c (which mimic the standard model particles as mentioned in Ftn. [T]), but these corrections are 
suppressed by the large masses of the a\, scalars (in the simple model, all masses are taken to be equal to 
M). Still, radiative corrections may provide a valuable alternative to direct searches if all masses of the 
new tJ-type particles are of the order of several TeV or more. 



5 



with e = ey for the scalar theory f l2.2ap as it stands and an appropriate constant value go of 
the dynamic g-field (or an appropriate value /iq of the corresponding chemical potential /i; 
see Refs. 0- J ^ for further discussion). 

The really new input for model fl2.1al) is the following Ansatz for the effective 
gravitational-coupling parameter: 

1 1 r . ^ 1 f ?/2 for T > t}^J , 

Kriq) = 7. go + 7. (g - go) ^ T/t};J - 1 = (2.2e) 



go/2 for T < 



{+) 

c,K ; 



2 ' 2 

with the step function 

n I 1 for X > , , , 

= , , ^ ' (2.2f) 
10 for X < . 

For the model universe to be discussed in Sec. Hill having a temperature decreasing with 
cosmic time, Kx is a nontrivial function of g above a critical temperature T^~^ and a 
constant below Tj~^. A possible physical realization of the corresponding "first-order phase 
transition" will be given in Sec. Ill CI Here, two brief remarks suffice. First, there may be 
hysteresis-type effects, and the suffix '(+)' on the critical temperature is to indicate that the 
transition is approached from the high-temperature side. Second, there are the following 
assumptions on the critical temperatures entering (I2.2cp and (I2.2ep : 



T,,, = 0(i?ew), T,;+) = 0(i?ew) , T,,g>T};^. (2.2g) 

The constant go in f l2.2ep . relevant at zero temperature, is proportional to the inverse of 
Newton's gravitational constant, specifically 

go = l/(87r Gn) = (^Pianck)' ~ (2.44 x 10^^ GeV)^ . (2.3a) 

Given the energy scale from and (^^, a single dimensionless parameter char- 



acterizes the theory, namely, the ratio of the two energy-density scales, 

^ = {Epii^^ck/ EewY , (2.3b) 
which is approximately 10^° for E^w ~ 2.44 TeV. 

C. Additional remarks 

The previous subsection has given the detailed description of the field-theoretic model 
(12. lap . Still outstanding is the promised physical realization of the gravitational-coupling 
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Ansatz fl2.2ep . The particular realization relies on the possibility of having symmetry restora- 
tion at low temperatures and symmetry breaking at high temperatures, the opposite of what 



is the case in most systems. This possibility has been discussed in, e.g., Ref. 20| and it is 
the easiest to just follow Example 3 of Sec. IV of that article. 

The argument, then, pro ceeds in four steps. First, start from a scalar theory with global 
0{n) X 0{n) symmetry [20|, where the scalar fields are denoted as xa and respectively, 
with both indices A and a running over 1, . . . , n [again, an inner-product notation will be 
used, X ■ X = J2a XaXa and r] ■ r] = J2a VaVa]- Second, take the parameters in the zero- 
temperature potential in such a way as to give the following symmetry-breaking pattern in 
a finite-temperature context: 

0{n) X 0{n) ^ 0{n) x 0{n - 1) , (2.4a) 

T=0 T=oo 

which implies that the rj scalars develop an expectation value at high temperatures. As 
indicated by the superscript on the arrow in fl2.4ap . the finite-temperature phase transition 
is arranged to be first-order.'^ Third, consider the following hypothetical interaction term in 
the effective Lagrange density: 

A-„fi+(,/2-A-„) ,^j;;;,;;,^„^ fi. (2.4b) 

where the fraction is close to 1 for field values x ' X '^''^ V ' V very much larger than (i^ew)^, 
which is the case relevant to temperatures far above the phase transition. Fourth, with a 
single ultraviolet parameter Kq = go/2, the finite-temperature behavior f l2.4ap for the term 
(I2.4bl) essentially gives the previous Ansatz (\2.2e\\ . Again, the phase-transition explanation 
of the postulated behavior fl2.2ep does not need to be taken literally: the origin of Kxiq) 
may very well have an entirely different origin, assuming the model to be relevant at all. 

It is clear that Eqs. fl2.ip -( l2l2l) . as they stand, only provide a phenomenological model. 
The main ingredient is the discontinuous phase-transition- type behavior of (I2.2ep . The 
model is simpler than the one used in Ref. 0] and, more importantly, entirely within the 
framework of g-theory (which provides a possible solution of CCPl). 



^ Specifically, the potential P(x, rj) is given on p. 3367, left column of Ref. 20| with mass-square parameters 
now positive and of order (i?cw)^, to which are added the following "cubic" terms: gi E^^ (x ■ xf^"^ and 
92 Ecw iv ■ ?/)^^^. The coupling constants are taken as in the fourth unnumbered equation on p. 3368, left 
column of the same reference, together with, for example, gi < and 52 < 0. These cubic terms and the 
finite-temperature corrections of order Tr]'^ can give a first-order (discontinuous) phase transition 19|. 



The first-order nature of the phase transition (|2.4ap may also have other origins, see Endnote [24] of 
Ref. 0- 
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III. COSMOLOGY 



A. Setup 



A spatially flat, homogeneous, and isotropic universe with scale factor a{t) will be consid- 



ered 
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For convenience, this cosmology will be called a Friedmann-Robertson- Walker 
(FRW) universe, even though, as will become clear in Sec. IIII CI the Friedmann equation is 
slightly modified when the electroweak kick sets in. 

The homogeneous matter content of this model universe consists of two perfect fluids 
(called type 1 and type 2), with energy density Pmi(^) from massive a scalars with an effective 
number of degrees of freedom A^i and energy density pm2 (t) from massless ip scalars with A^2 
degrees of freedom (see the first paragraph of Sec. Ill Al and the Appendix of Ref. 0] for the 
physical motivation of having Ni = N2 = 10^). In thermal equilibrium and without energy 
exchange, the type-2 energy density is given by pm2 = {N2 n / 30) T'^ . The temperature of 
the Universe can, therefore, be identified as approximately (pm2)^^^- 

For pmi = and pm2 ~ (-E'ew)^, the expansion of the Universe is governed by a Friedmann- 
type equation (see below) with a timescale set by 

tew = -E'pianck/ (-E'cw)^ , (3.1) 

in terms of the reduced Planck energy from ( I2.3ap . A value i?ew ~ Te V gives 1/t ew ~ meV. 



B. Frozen-electroweak-kick mechanism 



With the field-theoretic model of Sec. [TTl for an energy scale E^.^ ~ TeV, the basic steps 
of the frozen-electroweak-kick mechanism in a fiat FRW universe are as follows: 

(i) start from a standard radiation-dominated FRW universe at an ultrahigh temperature 
T with pv = Pmi = and pm2 ~ T'^ from the massless scalars ip (the ip scalars 
have standard electroweak interactions and are in thermal equilibrium, whereas the 
a scalars may have nonstandard interactions and are assumed to be initially absent, 
Pmi = 0); 

(ii) as the temperature T drops below ^ (with T^. ^ > -Ecw), the a"^ coupling of the 
scalar theory fl2.2ap generates a nonzero density of massive a scalars at cosmic times 
t around tc,g (with < ^ew); 

(iii) the presence of massive scalars a modifies the Hubble expansion rate H{t) = 
a{t)~^ da{t) / dt at t ~ tew; 
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(iv) the modified Hubble expansion rate kicks pv away from zero {gI, witli pv{t) ~ H{t)^ ~ 
(l/tew)"^ ~ (meV)"^ at t ~ 

(v) a nonzero value of pv remains when Kt from fl2.2el) is frozen to the constant value 
go/2 at a temperature T = Tc^k (with Tc^k ^ -^ew) or cosmic time t = tc,K (with 

tc, K ^ '^ew) ) 

(vi) the subsequent evolution is that of a standard A-FRW universe (here, with relativistic 
scalars ip, as the massive scalars a ultimately disappear).^ 

The phenomenological model of Sec. [Ill is, most likely, over-simplified, but may provide 
a benchmark calculation for a dynamically generated cosmological constant. Expanding on 
items (iv)-(vi) above, note that, at cosmic times t ~ tew, the frozen vacuum energy density 
Pv, remnant ~ (l/tew)"^ is negligible compared to the matter energy density pM2{tew) by a factor 
^ ~ 10^°. With Phnif) oc the tiny (but constant) vacuum energy density py, remnant 

only becomes dominant at very much later times, tyM-equai ~ 4w ~ 10^^ s, suggesting a 
possible solution of the so-called cosmic coincidence puzzle (see, e.g., Ref. jsl, [lol, 11|). 



In the next subsection, a preliminary numerical calculation is presented, which supports 
the scenario of the frozen-electroweak-kick mechanism outlined above. As this next subsec- 
tion is rather technical, it may be skipped in a first reading, except for a quick look at the 
numerical results in Fig. [1] 

C. Numerical solution 

With the timescale tew from (13. ip and the hierarchy parameter ^ from ( I2.3bll . it turns 
out to be useful to introduce the following dimensionless variables for the cosmic time, the 
Hubble axpa^ion .ate, the energy den^itie. ..d the , =Utt away from e,uihb„utn g 

r = (tew)"^i, h = te^H, (3.2a) 

rMn = {te^Y PMn , Ty = {te^Y PV = 3:^/2 , (3.2b) 

X = e (g/go - 1) , (3.2c) 

where n stands for the matter-species label [n = 1, 2), and the pv Ansatz f l2.2dp has been 
used. From now on, an overdot will denote differentiation with respect to r. 



^ Strictly speaking, such a standard A-FRW universe would not allow for further large contributions to the 
vacuum energy density at temperatures T < Tc^k ^ Ecw Naively, one expects a contribution of order 
(100 MeVY_ from quantum chromodynamics at T ^ 10^ MeV, but it has also been argued that this is not 
the case [231. 
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The relevant dimensionless ordinary differential equations (ODEs) for the model of Sec. |TT] 
are then given by 

{h + 2h'^)^x'^/2 + ^{rMi + rM2-3h^) -Sh'^xe -hxx = 0, (3.3a) 
r-Mi + (4 - kmi) hvMi = A21 rM2 - A12 tmi, (3.3b) 
fM2 + 4:hrM2 = -A21 rM2 + A12 Tmi, (3.3c) 
{3hx + 3h^x)e- xy2 + ^{rMi + rM2-3h^) =0, (3.3d) 
with the EOS function kmi{t) defined in Sec. A2 of Ref. p] and the effective step function 
e{r)=e[rM2{r)-r,,K], (3.4a) 



using definition (I2.2fl) . Here, the energy density rM2 of the massless scalars monitors 
the ambient temperature of the model universe and determines the moment when the g- 
dependence of K{q) changes from q/2 to go/2, as given by Ansatz fl2.2el) . For the benefit of 
the reader, the four ODEs in (13. 3p trace back to Eqs. (4.1) and (4.2a) of Ref. [sj evaluated for 
the Rnp-Ansatz fl2.2el) and Eqs. (A5b) and (A5d) of Ref. [7j adapted to the case considered. 

As discussed in Sec. Ill Bt the scalar interactions turn on below a certain critical tempera- 
ture, which corresponds to rM2 < fc,g in the cosmological context. The postulated behavior 
f l2.2cl) suggests the following coupling parameters in the ODEs (13.31) : 

2 

(3.4b) 



Ai2(r) = X9[rc,g -rM2] i 1 - \jrM2/rc,g 



A2l(l 



Ai2(r) exp 



TiN2 V'^ a(r) M 

30rM2(T-mm)/ ^(^min) ^ew 



(3.4c) 



with A oc {qoY- The argument of the exponential in f l3.4cp equals the negative inverse of 
the T/M expression (A3d) from Ref. 0], and Tmin is an arbitrary reference time before Tbcs 
to be introduced below. The exponential factor of (I3.4cp ensures that the ODEs (]3.3bp and 
f l3.3cp for Minkowski spacetime {H = 0) at a fixed temperature T < T^^g (rjv/2 < ^c^g) give 
an equilibrium ratio rMi/rM2 = exp[— M/T]. 

The new physics from Eqs. fl2.2cp and (I2.2ep is assumed to operate in a temperature 
window set by rM2(T) values between r^g and r^i^, with 



c,g 



(3.5) 



according to (|2.2gp. In this regime, the square bracket in fl3.3dp corresponds to the standard 



Friedmann equation {H^ oc ptot), to which are added two terms (proportional to q and q — qo) 
tracing back to the g-dependence of the gravitational-coupling parameter Q . 
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The ODEs (13.31) have two interesting analytic solutions. The first corresponds to a 
standard radiation-dominated FRW universe at high enough temperatures [with Ai2(t) = 
A2i(r) = 0], 



hir) = (1/2) , x{t) = rMi(r) = , rM2{r) = (3/4) > r,,, , 



(3.6a) 



and the second to a standard de-Sitter universe with constant vacuum energy density and 
without matter, 



h =^ X /6, h = x = 0, rMi = rM2 = 0. 



(3.6b) 



The numerical solution to be presented shortly will be seen to interpolate between these two 
analytic solutions. 

The hypothetical TeV-scale physics has a very large hierarchy parameter ^ ~ 10^*^ from 
(l2.3bD . and four ^ = oo equations turn out to be relevant for the phase of the electroweak 
kick [7.]. Specifically, there are two differential equations. 



3 (^h + 2h^) - -Kmi {3h^-rM2) 



{e)hx 

rM2 + 4 /i rM2 + A21 rM2 = A12 (3 - rM2) , 



{i-e)h 



,x <x 



and two algebraic equations. 



X 



- kmi (3 



f M2 



rM2) 



+ (1 



f M2) 



rM2=rc, K 



(3.7a) 
(3.7b) 

(3.8a) 
(3.8b) 



with 6 defined by (13. 4p . Recall that kmi is a functional of /i(r) = a(r)^^ da{T)/dT, explicitly 
given by Eqs. (A3a)-(A3d) in Ref. [7|. For the numerical analysis of Eqs. fl3.7p - (l3.8l) . the 
relevant parts of the equations for the 6 = 1 (high-temperature) phase have been indicated 
by square brackets, and those for the = (low-temperature) phase by curly brackets, 
keeping the two generally valid equations without such brackets. 

Returning to general values of ^, the numerical solution of the ODEs (13.30 has been 
obtained for the case with Ni = N2 = 10^ and an equal mass M of all type-1 particles 
(i.e., the case-B mass spectrum in the terminology of Sec. A2 of Ref. [t']). As mentioned 
above, at temperatures above the critical temperature T^^g, it is possible to have a standard 
radiation-dominated FRW universe (I3.6al) with all type-2 particles in thermal equilibrium 
(the type-1 particles are assumed to be absent in the early phase, see Sec. IIII B[) . Hence, 
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FIG. 1: Numerical solution of the dimensionless ODEs (j3.3p with EOS function kmi{t) defined 



in Sec. A2 of Ref. 



a 



and further parameters (j3.4p . The panels are organized as follows: the four 
basic dynamic variables [/i(t), rMi^T), rM2{T), and x{t) as defined by ()3.2p ] are shown on the top 
row and secondary or derived quantities on the bottom row. The dashes lines in the panels of the 
third column correspond to rM2 and A21/A. The main result is the nonzero remnant value 
of the dimensionless gravitating vacuum energy density ry = shown in the top right panel. 

The model parameters are A, Tc^g, Tc^k} = {lO'^, 10^, 12, 3}. The ODEs are solved over the 
interval [rmm-, ''"max] = [0.01, 10] with the following boundary conditions (j3.9p at r = rtcs = 0.25 
corresponding to rM2(Tbcs) = ^'c.g : {x, h, a, vmi, rM2} = {0, 2, 1, 0, 12}. 

the appropriate boundary conditions on the four dynamical variables at a time r = Tbcs are 
/i(i"bcs) = ^(rbcs)^\ (3.9a) 
rMiir^cs) = 0, (3.9b) 
'"M2(7-bcs) = 3 [/^(rbcs)]^ (3.9c) 
x(rbcs) = 0. (3.9d) 

The precise value of Tbcs is irrelevant as long as it is sufficiently small, with rM2{Thcs) > i^c,g 
[physically interpreted as T(tbcs) > ^c,g]- Furthermore, choose for A the value 10^, making 
rMiir) decrease significantly before the Kt transition at rM2 = r^K is reached (the new 
physics then operates in a relatively narrow temperature interval near T ~ -Eow)- Other 
values A > 10^ give similar numerical solutions. 

The numerical results are shown in Fig. [1] The ry panel, in particular, shows the narrow 
new-physics window with the critical temperature fl2.2cp at r = 0.25 (from rM2 = ^c,g) and 
the freezing of the gravitational-coupling parameter fl2.2ep at r ~ 0.45 (from rM2 = ^c,k)- 
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TABLE I: Asymptotic values of the dimensionless gravitating vacuum energy density ry(r) for 
various hierarchy parameters ^ ranging from 10^ to 10^. All other parameters are given in the 
caption of Fig. [TJ The entry for ^ = 00 has been calculated from Eqs. (j3.7p and (j3.8p . The 
numerical accuracy is estimated at ±1 in the last digit shown. 





10^ X rv{oo) 


10^ 


7.310 


10^ 


5.380 


10^ 


2.028 


10^ 


2.376 


108 


2.376 


00 


2.376 



The calculated numerical value of ry^ remnant is approximately 2.4 x 10^^. For different values 
of Tc^K than chosen in Fig. [H while keeping the other parameters the same, the numerical 
values of ry^ remnant wiU, of coursc, bc Icss than the maximal value of ry(r) shown in the top 
right panel of the figure, that is, ry_ remnant ^ 1-5 x 10^^. Incidentally, the peak of the rv'(r) 
curve at r = Tpeak ~ 0.275 looks sharp in the plot but is really a concave parabola with a 
large negative second derivative, ry(rpeak) ~ —10^. 

Table H] presents the numerical values for ry(oo) = 1101^-^00 ^vi'^) for various hierarchy 
parameters ^. The ry(r) solution for ^ < 10^ has, in fact, significant oscillations superposed 
on the smooth curve shown in the bottom right panel of Fig. [T], which explains the somewhat 
erratic behavior of the first three entries in Table [H Based on the analysis of Ref. 0], the 
results for C, = 00 can be expected to give a close approximation to those for C, ~ 10^°, which 
corresponds to the physically relevant case according to (I2.3bl) . 



D. Analytic result for ^ = 00 



The ^ = 00 equations (13. 7p and (13.80 immediately give an analytic expression for the 
asymptotic value of the dimensionless vacuum energy density. 



lim ry(r) 



^=00 



/tMllTfrceze, 



[3 hi 



Tfr, 



^M2(Tfrcozc) = ''c, K 



(3.10) 



where the EOS function kjv/i(t) is determined by the dimensionless Hubble expansion rate 
/i(r) and VMni'^) stands for the dimensionless energy density of matter component n [all 
dimensionless variables are defined in (13. 2p ]. 
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Expression fl3.10p is formal, because the numerical solution of the two nonlinear ODEs 
fl3.7ap and (13.7b|) for = 1 is needed to determine the value of, for example, /i(rfi.eeze)- But 



the analytic expression (13.101) does clarify the essential physics involved: the EOS function 
KMiit) multiplied by the corresponding energy density of massive particles {rjui = 3 h?—rM2) 
and the freezing of K{t) = q{t)/2 to the constant value go/2 at a cosmic temperature of the 
order of E^^ ~ TeV. 

IV. CONCLUSION 



I2.1l) - fl2.2p has been presented, which re- 



in this article, a simple field-theoretic mode 
mains entirely within the framework of g-theory 2- J] and does not require unnaturally small 
coupling constants (the energy scales £'ew ~ TeV and -Epianck ~ 10^^ TeV are considered to 
be given j^). As summarized in Sec. IIII Bl the model generates via the electroweak-kick 
mechanism |6| an effective cosmological constant Ace (remnant vacuum energy density pv), 
which is consistent with the value A"*^*^ ~ 1 x 10~^^ gcm^'^ ~ (2 x 10^'^ eV)"^ from obser- 
vational cosmology [ill, [21|. In addition, having Aes ~ (-E'ew)^/(-£'Pianck)'^ provides a natural 



explanation [5| of the fact that the orders of magnitude of the energy densities of vacuum, 
matter, and radiation are approximately the same in the present Universe (also known as 
the triple cosmic coincidence puzzle). 

With the calculated dimensionless vacuum energy density ry(oo) ~ 2.4 x 10^^ from 
TablelH the required energy scale E^^ of the new physics is approximately 4.7 TeV, according 
to Eq. (5.2) of Ref. 0]. However, the main focus of the present article is not on numerical 
estimates (plenty have been given in Ref. Q), but rather on the physical content of a theory 
capable of generating the observed cosmological "constant" of our Universe. 

In that spirit, the most interesting result of this article is the observation that the pro- 
posed model involves one crucial ingredient, namely, the discontinuous phase-transition- type 
behavior of the dependence of the gravitational coupling f l2.2ep on the quinta-essentia field 
q. (The need for some form of singular behavior in order to freeze pv[q{t)] has also been 
emphasized in the second remark of Sec. Ill E of Ref. [7].) The main task is, therefore, 
to find the rationale for this phase-transition type of behavior or for a different effect with 
the same result of freezing part of the vacuum energy density below a certain cosmic tem- 
perature scale. Furthermore, the simple model (or a suitable generalization of it) needs to 
be embedded in the complete solution of the cosmological constant problem, which is still 
outstanding. 
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